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1. Introduction 

The renormalization group (RG) based technique of "lattice matching" via block transforma- 
tions relates physical quantities on different lattices. It thus provides a method for computation of 
a physical quantity at different lattice spacings (couplings). To apply it one needs to implement RG 
block transformations on the approach to the Wilsonian Renormalized Trajectory (RT). This can be 
done in various ways. One way is numerical implementation of the RG blocking by Monte Carlo 
RG (MCRG) techniques. This is the method that has mostly been used in the literature. Another 
approach is to implement blocking by explicit RG recursion relations that can, to varying degree, 
be explicitly carried out by numerical-analytical means. This is the method followed here. Specif- 
ically, we employ explicit RG recursion relations of the "potential moving" type. These block 
transformations (decimations) are of course approximate but can, in principle, be systematically 
improved. They turn out to be surprisingly effective. In the following we apply lattice matching 
of these decimations to obtain critical couplings and string tensions for the SU (2) and SU (3) pure 
gauge theories. A more detailed account has appeared in [JTJ]. 

2. RG blocking recursions and lattice matching 

RG blocking recursions We start by assuming a general plaquette action A p (U p ,n) on lattice 
of spacing b n a given in terms of the character expansion of its exponential: 

exp {-A p (U p ,n)) = £ djFj(n)Xj(U p ) . (2.1) 



The sum is over all inequivalent irreducible representations labeled by j, with characters Xj °f 
dimension dj. The action itself is, of course, completely specified by the set of Fj(n) coefficients, 
and vice versa, and of the general form: 



A p (U p ,n) =£ jH n )^[X0 P )+Xj{U p 1 )] 

j 



(2.2) 



with lj = 1 for self-conjugate and lj = 2 for non-self-conjugate representations. (For SU{2), in 
particular, lj = 1 for all j.) 

It is useful to define an effective coupling gto characterizing a given action of the form (2.2): 
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Here {?} are the SU(N) generators and m a unit vector. ((2.3) is of course independent of the 
direction m). In the perturbative regime this reduces to the usual definition of gauge coupling. In 
the non-perturbative regime any definition of a 'coupling' is of course some scheme-dependent 



choice. We adopt ( |2.3| ) to track the RG recursion flows; it provides an efficient parametrization of 
the renormalized trajectory below. 

The lattice block step b" a — > b" +1 a may now be formulated as a prescription for the character 
expansion coefficients Fj(n + 1) in terms of the Fy(n)'s: 
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(2.4) 
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To complete the prescription we must specify the renormalization parameters £, r. We take 

f = b\l-cg^ 2 ] (2.5) 
r = b \\-cg (n)2 \ (2.6) 

with c an adjustable decimation parameter to be tuned for optimization as explained below. It is 
convenient to work with normalized coefficients fj = Fj/Fq by factoring out the trivial representa- 



tion coefficient in (2.1 ). Effective couplings (2.3) are also conveniently computed directly in terms 
of the {fj}. 

RG flows and lattice matching Consider now a general lattice system described by an action 
A(K) with set of couplings K = {Kj}. Successive RG blockings by a scale factor b generate a flow 
in action space: K — > 

i) £-(2) >. #(») . . . f w here = {K^} denotes the couplings 

after n blocking steps. If flows from K and K' reach the same point on RT after n and n' steps, then 




K K 



{Kj} 



Figure 1: Flows towards the RT from two different starting couplings K and K'. 
the corresponding lattice con-elation lengths and spacings a, a' are related as 

= a' = Z?("-"V (2.7) 

To identify such pairs of couplings we need ascertain that after n and n' RG steps, respectively, 
the same point is reached on the RT. This can be done in either of two ways: (i) show that the cor- 
responding actions coincide: A(K^) =A(K'^). This requires that one obtain the blocked action 
at each step; or (ii) show that the expectations of every operator, measured after performing the 
corresponding number of blocking steps from the initial two actions, agree. Either way, identifying 
such pairs (K, n), and (K', n') is referred to as two-lattice matching 

If blockings are performed numerically by MCRG, the second method appears easier to use. 
Obtaining the blocked action can be difficult, whereas it is possible, at least in principle, to generate 
a Boltzmann-weighted configuration ensemble for the blocked action by instead blocking the con- 
figurations of an ensemble generated from the original action. These can then be used to measure 
observables In practice, of course, due to lattice size limitations, only a rather small number 
of block steps is possible by MCRG, so getting close enough to the RT is not guaranteed. As a 
general observation, the location of the fixed point being block definition dependent, appropriate 
fine-tuning of any decimation free parameters can be crucial for achieving rapid approach in few 
steps. 



3 



Critical couplings from two-lattice matching ofRG decimations 



E. T. Tomboulis 



Here we employ two-lattice matching with RG block transformations implemented by the 
recursions (2.4) described above. They can be explicitly evaluated to any desired accuracy on 
lattices of any size, so no inherent limitations due to finite size arise. The blocked action resulting 
after each RG step is explicitly obtained, so it can used to ascertain approach to the RT and perform 
two-lattice matching. The transformations contain one parameter (cf. (2.5) - (2J5)), which should 
be fixed for matching optimization. 

A basic feature of our decimations is that, regardless of the choice of the initial plaquette 
action, a single step suffices to generate an action of the form (2.2) generally containing a large 
(infinite) set of representations. This is important as flow in such a large-dimensional interaction 
space makes it possible to avoid getting stuck at (finite-dimensional) lattice artifact boundaries. 
Furthermore, MCRG construction of blocked actions [Q] shows that one-plaquette terms with a 
large number of characters are the most relevant action terms for long-scale dynamics. This is 
precisely the type of action resulting from our decimations, and may be the reason for their apparent 
efficacy in computing long-distance dynamics observables as seen below. 

In the following the starting action for our decimations (n = 0) will always be taken to be the 
fundamental representation Wilson action. (Other choices such as mixed actions containing several 
representations can be treated in exactly the same way.) One finds that the flow under successive 
decimations reaches a unique RT irrespective of such a choice, though of course the number of 
steps needed to reach it depends on the initial point in action space. With the fundamental Wilson 
action as the starting action the approach is found to be very rapid as illustrated in Fig. 0. 




Figure 2: Flow from the SU (2) fund. Wilson action with j3 = 4 (green dots) and j3 = 2.5 (red dots) showing 
rapid approach to the RT. First three non-trivial (normalized) expansion coefficients shown. 



The effective coupling (2.3) provides a good way to label points along the RT. If, starting from 
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some Wilson action coupling /3, after n steps the point j3 W (j3) lies on the RT, subsequent RG steps 
generate a sequence of points j3'" +1 \ j8^" +2 \ • • • hopping along the RT. With scale factor b = 2, and 
for all large and intermediate values of j3 W, the effective beta function is varying slowly enough for 
a linear interpolation to provide an excellent approximation to the RT points lying between pairs of 
neighboring points fi( n \ fi( n+l \ 

The decimations become exact in the j8 — > °° limit. Computation of the step scaling function 
(beta function) from the decimations in the weak coupling scaling region reproduces the perturba- 
tion theory prediction to within 2% — 3%. We next use them to obtain critical couplings and string 
tensions by lattice matching for the SU (2) and SU (3) gauge theories. 

3. Critical couplings and string couplings by two-lattice matching 
3.1 Critical couplings 

At physical temperature T = 1 /aN t , lattice with temporal extent N T , spacing a and lattice with 
N' T , a' are related by: 

a = —a. (3.1) 
If after blocking n and n' times, respectively, the two flows reach the same point on the RT, using 



(2.7) this implies 



n-rt = 1og b (j£\ . (3.2) 
So, at T = T c one has 

pto(Pc(N*)) = P W) (Pc(K))- (3-3) 
This suggests the following simple matching procedure. 

(i) Assuming j3 f (Af T ) known for one N x , take 

n = \og b N z + m 

ri =\og b N[+m (3.4) 

Integer m = 0, 1, . . . is chosen so n,n' large enough to be on the RT. (If the so-chosen n or/and n' 
turn out to be non-integer, one performs [n] and [n] + 1 steps, where [n] is the nearest integer to the 
chosen n from below, and uses interpolation for the RT points in-between as mentioned above.) 

(ii) With n, n' and j3 c (A^) given, solve ( pT3| ) for f} c (N' T ). This means that the starting point of the 
flow on the N' T lattice is adjusted to satisfy (^3|). 



3.2 String tensions (T = 0) 

A similar procedure allows one to obtain string tensions by matching. Assume that two RG 

flows from starting Wilson action at j8o and fii end up at the same RT point after no and n\ steps, 
respectively. Then one has 

pM(Po) = pW(p l ) (3.5) 

and 

a l ^fo = b {n °-' n) ao^fo (3.6) 



Suppose aoy/o known. Choose no large enough to be on the RT. Determine n\ so that (3.5) is 



satisfied. a\y/a is then obtained directly from (3.6) 
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3.3 Numerical results 



It is important to maintain high accuracy in working with the expansions (2.1) under blocking 
iteration. For SU(2) we typically use fifty group characters in the expansions (2.1). This implies 
for, say, /3 = 5 omitted higher character coefficients fj = Fj/Fq, and accompanying bounds on the 
series remainder, of the order of 10~ 45 . For SU (3) we truncate (2d) at characters j = (p,q) with 
p>20,q> 20; this implies remainders at /3 = 10 of less than 10~ 12 . 

The scale factor is always taken to be b = 2. The adjustable parameter in the decimation 
recursions ( p74j ) - (2.6) is c, which we set at c = 0.10 in the case of SU (2) and c = 0.24 in the case 
of SU (3). With no other parameters present, straightforward numerical evaluation of the recursion 
relations can then be carried out. 

We take one value of p c (N t ) from MC data, which serves to fix the scale and apply the pro- 
cedure above to obtain critical coupling values for other lattices. Results for SU (2) are shown in 
Table [l|. Two sets of computed f3 c values are shown in Table [j] (columns 1 and 2) corresponding to 
two different choices of the MC data point (underlined entries). The table also shows comparison 
with the values obtained by MC simulation [||] - [Q] in each case (column 3). The agreement is 
very good - at the 2% — 3% level. Results for critical couplings in the SU (3) gauge theory are 
displayed in Table ^[ Agreement with MC simulation data []|] is again very good, typically within 
a few percent. 





ft 


ft 


j3 c (MC) 


3 


2.1875 


2.1957 


2.1768(30) 


4 


2.2909 


2.2991 


2.2991(02) 


5 


2.3600 


2.3683 


2.3726(45) 


6 


2.4175 


2.4258 


2.4265(30) 


8 


2.5097 


2.5180 


2.5104(02) 


12 


2.6355 


2.6440 


2.6355(10) 


16 


2.7275 


2.7361 


2.7310(20) 


32 


2.9487 


2.9574 





Table 1: Critical couplings j5 c (N T ) for SU(2) computed from lattice matching of decimations. Column 1 
and 2 show the values obtained for two different choices (underlined entries) of the one data point taken 
from MC data (see text). Column 3 shows the values from MC simulations for comparison. 



String tensions in SU (3) obtained by the method above are displayed in the same format in 
Table [|. Very similar results are obtained for SU (2) Good agreement with MC data - [§] is 
again obtained in all cases. 

Consideration of fermionic observables by similar RG recursion methods is a rather more 



demanding proposition. Some preliminary attempts are reported in fllOfl. 

This work was partially supported by the NSF under NSF-PHY-0852438. 
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